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Abstract 

We use Rice's formulas in order to compute the moments of some 
level functionals which are linked to problems in oceanography and op- 
tics. For instance, we consider the number of specular points in one or 
two dimensions, the number of twinkles, the distribution of normal angle 
of level curves and the number or the length of dislocations in random 
wavefronts. We compute expectations and in some cases, also second mo- 
ments of such functionals. Moments of order greater than one are more 
involved, but one needs them whenever one wants to perform statistical 
inference on some parameters in the model or to test the model itself. 
In some cases we are able to use these computations to obtain a Central 
Limit Theorem. 
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1 Introduction 

Many problems in applied mathematics require to estimate the number of points, 
the length, the volume and so on, of the level sets of a random hmction W{x), 
where x G R'', so that one needs to compute the value of certain functionals of 
the probability distribution of the size of the random set 

Cj{xi,uj) := {x e A : T4^(x,t^) = u}, 

for some given u. 



Let us mention some examples which illustrate this general situation: 

*Universite de Toulouse, IMT, LSP, F31062 Toulouse Cedex 9, France. Email: azais@cict.fr 
^Escuela de Matematica. Facultad de Ciencias. Universidad Central de Venezuela. A. P. 

47197, Los Chaguaramos, Caracas 1041-A, Venezuela. Email: jose.leon@ciens.ucv.ve 

tCentro de Matematica. Facultad de Ciencias. Universidad de la Repiiblica. Calle Igua 

4225. 11400. Montevideo. Uruguay, wschebor@cmat.edu. uy 



1 



• The number of times that a random process {X{t) : t G M} crosses the 
level u: 

Generally speaking, the probability distribution of the random variable 
(u) is unknown, even for the simplest models of the underlying process. 
However, there exist some formulas to compute E(iV^) and also higher 
order moments. 

• A particular case is the number of specular points of a random curve or a 
random surface. 

Consider first the case of a random curve. We take cartesian coordinates 
Oxz in the plane. A light source placed at {0,hi) emits a ray that is 
reflected at the point {x,W{x)) of the curve and the reflected ray is reg- 
istered by an observer placed at (0, /12). 

Using the equality between the angles of incidence and reflexion with re- 
spect to the normal vector to the curve - i.e. N{x) — {—W'{x),l) - an 
elementary computation gives: 

W'ix) = ^^HTIZ^ (1) 

x(r2 - n) 



where := — W{x) and := y^x^ + a^, i=l,2. 



The points (x, W{x)) of the curve such that x is a solution of ([T]) are called 
"specular points". We denote by SPi{A) the number of specular points 
such that X G A, for each Borel subset A of the real line. One of our aims 
in this paper is to study the probability distribution of SPi{A). 

• The following approximation, which turns out to be very accurate in prac- 
tice for ocean waves, was introduced long ago by Longuet-Higgins (see [13] 
and [14]): 



Suppose that hi and /i2 are big with respect to W{x) and x, then 
ai + x^/{2ai) + 0{h~^). Then, (jT]) can be approximated by 



W'ix)^^^^±^^^^±^^kx, (2) 
2 aia2 2 n.i«2 



where 



k 



— -(— — 



Denote Y{x) :— W'{x) — kx and SP2{A) the number of roots of Y{x) 
belonging to the set A, an approximation of SPi {A) under this asymptotic. 
The first part of Section [3] below will be devoted to obtain some results 
on the distribution of the random variable 5*^2 (K) ■ 
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• Consider now the same problem as above, but adding a time variable t, 
that is, W becomes a random function parameterized by the pair {x,t). 
We denote Wx,Wt, W^t, ■■■ the partial derivatives of W. 

We use the Longuet-Higgins approximation ([2]), so that the approximate 
specular points at time t are (x, W{x, t)) where 

Wxix, t) = kx. 

Generally speaking, this equation defines a finite number of points which 
move with time. The implicit function theorem, when it can be applied, 
shows that the x-coordinate of a specular point moves at speed 

dx _ Wxt 
dt Wxx - k ' 

The right-hand side diverges whenever W^x — k — 0, in which case a flash 
appears and the point is called a "twinkle". We are interested in the 
(random) number of flashes lying in a set A of space and in an interval 
[0, T] of time. If we put: 

then, the number of twinkles is: 

TWiA, T) := ^{{x, <) e A X [0, T] : Y{x, t) = 0} 

• IjCtW : Q dW^ -> W^' with d> d' he a random field and let us define 
the level set 

C^(u) = {xeQ:W^(x)-u}. 

Under certain general conditions this set is a (d — d')-dimensional manifold 
but in any case, its {d — (i')-dimensional Hausdorff measure is well defined. 
We denote this measure by Ud-d'- Our interest will be to compute the 
mean of the crrf-d'-nieasure of this level set i.e. Wip-d-d'iCq {vl))] as well 
as its higher moments. It will be also of interest to compute: 

E[/ Y{s)dad-d'{.s)]. 

where Y{s) is some random field defined on the level set. Cabafia [7], 
Wschebor [19] {d' = 1) Azais and Wschebor [4] and, in a weak form, 
Zahle [5D] have studied these types of formulas. See Theorems [S] and [S] 

• Another interesting problem is the study of phase singularities, disloca- 
tions of random wavefronts. They correspond to lines of darkness, in light 
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propagation, or threads of silence in sound [B]. In a mathematical frame- 
work they can be define as the loci of points where the amplitude of waves 
vanishes. If we represent the wave as 

VF(x,i) = ^(x,i) +ir;(x,t), where x e M"* 

where ^, rj are independent homogenous Gaussian random fields the dislo- 
cations are the intersection of the two random surfaces ^(x, t) — 0, ri{x, t) = 
0. We consider a fixed time, for instance t = 0. In the case d — 2 we will 
study the expectation of the following random variable 

#{xe^: e(x,0) -r7(x,0) =0}. 

In the case d = 3 one important quantity is the length of the level curve 

C{^eS: e(x,0) =77(x,0) =0}. 

All these situations are related to integral geometry. For a general treatment 
of the basic theory, the classical reference is Federer's "Geometric Measure The- 
ory" 0. 

The aims of this paper are: 1) to re-formulate some known results in a 
modern language or in the standard form of probability theory; 2) to prove 
new results, such as computations in the exact models, variance computations 
in cases in which only first moments have been known, thus improving the 
statistical methods and 3) in some case, obtain Central Limit Theorems. 

The structure of the paper is the following: In Section [5] we review without 
proofs some formulas for the moments of the relevant random variables. In 
Section [3] we study expectation, variance and asymptotic behavior of specular 
points. Section[3]is devoted to the study of the distribution of the normal to the 
level curve. Section [5] presents three numerical applications. Finally, in Section 
[S] we study dislocations of wavefronts following a paper by Berry & Dennis [B] . 



Some additional notation and hypotheses 

Xd is Lebesgue measure in R'', ad'{B) the d'-dimensional Hausdorff measure of a 
Borel set B and M'^ the transpose of a matrix M. (const) is a positive constant 
whose value may change from one occurrence to another. 

If not otherwise stated, all random fields are assumed to be Gaussian and cen- 
tered. 

2 Rice formulas 

We give here a quick account of Rice formulas, which allow to express the 
expectation and the higher moments of the size of level sets of random fields 
by means of some integral formulas. The simplest case occurs when both the 
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dimension of the domain and the range are equal to 1, for which the first results 
date back to Rice [17] (see also Cramer and Leadbetter's book [8]). When 
the dimension of the domain and the range are equal but bigger than 1, the 
formula for the expectation is due to Adler [Tj for stationary random fields. For 
a general treatment of this subject, the interested reader is referred to the book 
[1], Chapters 3 and 6, where one can find proofs and details. 

Theorem 1 (Expectation of the number of crossings, d = d' = 1) LetW — 
{W{t) : t ^ 1} , I an interval in the real line, be a Gaussian process having C^- 
paths. Assume that Var(T/F(t)) 7^ for every t Cz I . 
Then: 

E{N^{u)) = jE{\W\t)\\Wit) = u)pw(t){u)dt. (4) 

Theorem 2 (Higher moments of the number of crossings, d — d' — l) Let 

m > 2 be an integer. Assume that W satisfies the hypotheses of Theorem\^and 
moreover, for any choice of pairwise different parameter values ti,...,tjn G / 
the joint distribution of the k-random vector {W{ti), ...,W{tjn)) has a density 
(which amounts to saying that its variance matrix is non- singular). Then: 

E(7Vf (u)(Arf (u) - l)...{N^{u) - m + 1)) 




E( l[ \W'{t,)\\W{ti) = ... = W{trn) = u)pwiu),...,Wit,„){u,-,u)dti...dt 



(5) 

Under certain conditions, the formulas in Theorems [T] and [2] can be extended 
to non-Gaussian processes. 



Theorem 3 (Expectation, d = d' > 1) Let W : A C M.'^ ^ M.'^ be a Gaussian 
random field, A an open set o/R"^, u a fixed point in M."^ . Assume that 

• the sample paths of W are continuously differentiable 

• for each t Cz A the distribution ofW{t) does not degenerate 

• P({3t e A : W{t) = u , det(T^'(t)) = 0}) = 
Then for every Borel set B included in A 

E(iV]f (u)) = / E[|dctW^'(t)||VK(t) =u]pv^(t)(u)dt. 
Jb 

If B is compact, both sides are finite. 

The next proposition provides sufficient conditions (which are mild) for the third 
hypothesis in the above theorem to be verified (see again [J, Proposition 6.5). 
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Proposition 1 Under the same conditions of the above theorem one has 
P({3t e A : W{t) = u , det(VK'(t)) = 0}) = 

if 

• Px{t){^) ^ for all X in some neighborhood ofu, 

• at least one of the two following conditions is satisfied 

a) the trajectories of W are twice continuously differentiable 
b) 

a{S)^ sup P{|detV7'(t)| < (5|T4^(t) = x} ^ 
xev{u) 

as S where V{u) is some neighborhood ofu. 

Theorem 4 (m-th factorial moment d = d' > 1) Let m > 2 be an integer. 
Assume the same hypotheses as in Theorem\^ except for (Hi) that is replaced by 

(Hi') for ti, ...,tm G A distinct values of the parameter, the distribution of 

{W{ti),...,Wit^)) 

does not degenerate in (R'')™. 

Then for every Borel set B contained in A, one has 
E [(7V^(u)) (iV^(u) - l)...(A^^(u) - m + 1)] 

= / E(\{\Aei{W\t,))\\W{ti) = ...^W{t„,)^u^ 

Piy(ti),...,H'(t„)(u, ...,u)fiti...dt,„, (6) 

where both sides may be infinite. 

When d > d' we have the following formula : 

Theorem 5 (Expectation of the geometric measure of the level set. d > d') 

Let W : A ^ be a Gaussian random field, A an open subset o/M'^, d > d' 
and u G M"^ a fixed point. Assume that: 

• Almost surely the function t ^ W(t) is of class . 

• For each t €z A, W{t) has a non-degenerate distribution. 

• P{3t e A, W{t) = u, W'{t) does not have full rank} = 
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Then, for every Borel set B contained in A, one has 



¥.{<Jd-d'{W,B))= /" E([det(W^'(t)(M^'(t)f )]'/'|M^(t) = u) Pwit){^)dt 



J B 



(7) 



// B is compact, both sides in are finite. 



The same kind of result holds true for integrals over the level set, as stated 
in the next theorem. 

Theorem 6 (Expected integral on the level set) Let W he a random field 
that verifies the hypotheses of Theorem\^ Assume that for each t £ A one has 
another random field : V ^ R", where V is some topological space, verifying 
the following conditions: 

• Y*{v) is a measurable function of {uj,t,v) and almost surely, {t,v) ~^ 

is continuous. 

• For each t £ A the random process (s, w) (W(s),y*(w)) defined on 
W X V is Gaussian. 

Moreover, assume that g : A x C(V,M") M is a bounded function, which is 
continuous when one puts on C(V,M") the topology of uniform convergence on 
compact sets. Then, for each compact subset B of A, one has 



3 Specular points and twinkles 
3.1 Number of roots 

Let W(t) : M'^ ^ K'' be a zero mean stationary Gaussian field. If W satisfies 
the conditions of Theorem [3] one has: 






(9) 



where 



^0 

H being the spectral mesure of W. 

Formula ([9]) is in fact the one S.O. Rice wrote in the 40 's see [17] . 
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3.2 Number of specular points 

We consider first the one-dimensional static case with the longuet-Higgins ap- 
proximation ^ for the number of specular points, that is: 

SP2{I) = e I ■■ Y{x) = W'{x) -kx^O} 

We assume that the Gaussian process {14^(a:) : a; e K} has paths and 
Var(VF'(a;)) is constant equal to, say, v"^ > 0. (This condition can always be 
obtained by means of an appropriate non-random time change, the "unit speed 
transformation") . Then Theorem 1 applies and 

E{SP2iI))^ [Ei\Y'ix)\\Yix)^0)pYi,)iO)dx^ U{\Y' {x)\)-^{ — )dx 

J I J I V V 

G{-k.a{x))-^i — )dv, (10) 

/ V V 

where a'^{x) is the variance of W"{x) and G(/i, cr) := E(|Z|), Z with distribution 

For the second equality in (|10|) , in which we have erased the condition in the 
conditional expectation, take into account that since Var(W^'(x)) is constant, 
for each x the random variables W'{x) and W"{x) are independent (differen- 
tiate under the expectation sign and use the basic properties of the Gaussian 
distribution) . 

An elementary computation gives: 

G(/i, a) = ^i[2^^l/a) - 1] + 2aM'^), 

where ip{.) and $(.) are respectively the density and the cumulative distribution 
functions of the standard Gaussian distribution. 

When the process W{x) is also stationary, — A2 and cr'^(x) is constant 
equal to A4. If we look at the total number of specular points over the whole 
line, we get 

E(5P2(R)) = ^%^ (11) 
k 

which is the result given by [13] (part II, formula (2.14) page 846). Note that 

this quantity is an increasing function of -^^-^ . 

Since in the longuet-Higgins approximation A: ~ 0, one can write a Taylor 
expansion having the form: 



Let us turn to the variance of the number of specular points, under some 
additional restrictions. First of all, we assume for this computation that the 
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given process {W{x) : x S R} is stationary with covariance function 
E{W{x)W{y)) = T{x — y). T is assumed to have enough regularity as to per- 
form the computations befow, the precise requirements on it being given in the 
statement of Theorem [T] 

Putting for short S SP2{"&), we have: 

Var(S') = E{S{S - 1)) + E(S') - [HS)f (13) 
The first term can be computed using Theorem [21 

E(5(5-l))=/ / ¥.{\W"{x)~k\\W'\y)-k\\W'{x)^kx,W'{y)^ky) 
■Pw'ix),w'{y){kx,ky) dxdy (14) 

where 

, , 1 r lk^{X2X^ + 2T"^{x-y)xy + X2y^h 
Pw^i.,^wiy,{kx,ky) - ,^^,._r-,(^_^) e^P [-2 Xl-r"Hx-y) \ 

(15) 

under the additional condition that the density (flS)) docs not degenerate for 
x^y. 

For the conditional expectation in (|14p we perform a Gaussian regression of 
W"{x) (resp. W"{y)) on the pair {W'{x), W'{y)). Putting z — x — y, we obtain: 

W"(x) = 9y{x) + ay{x)W'(x) + by{x)W'{y) 

r"'(z)r"(z) 



ay{x) = - 
by{x) = - 



A2-r"2(z) (16) 

A2r"'(z) 
xl-r"Hz)' 



where Oy{x) is Gaussian centered, independent of {W {x),W' (y)). The regres- 
sion of W"{y) is obtained by permuting x and y. 



The conditional expectation in (jl4p can now be rewritten as an unconditional 
expectation: 

EJ n ,,T./"^_Mi , ^"{z)x + X2y^ ^ ,^^,„, ^,r, , T"{-z)y + X2x 



;{|g,(x)-fcr"(z)[l+ ^,^"J^^+J^ ]||g.(^)-fcr"(-z)[l 



A2-r"2(z)J '"^ ' "'L- A2-r"2(z) 

(17) 

Notice that the singularity on the diagonal x = y is removable, since a Taylor 
expansion shows that for z w 0: 

r'"W[i.™^] = i^.(=.0(.')). (18) 

One can check that 

a'{z) = E((0,(x))2) = E((0.(y))2) = A4 - (19) 
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and 

E{9yix)9Ay)) = r(^)(z) + ^^^^S^. (20) 

Moreover, if Ag < +00, performing a Taylor expansion one can show that as 
z Ki one has 

(21) 

and it follows that the singularity at the diagonal of the integrand in the right- 
hand side of (11411 is also removable. 



We will make use of the following auxiliary statement that we state as a 
lemma for further reference. The proof requires some calculations, but is ele- 
mentary and we skip it. The value of H{p; 0, 0) below can be found for example 
in 0, p. 211-212. 

Lemma 1 Let 

H{p;t,,iy)=Ei\^ + t,\\rj + iy\) 
where the pair (^,77) is centered Gaussian, E(^^) ~ ^(77^) = 1, IE(C?/) — p. 



Then, 
where 

and 



H{p; ly) = H{p; 0, 0) + i?2(p; /i, 1^) 
H{p- 0,0)^-^/1- p^ + — arctan ^ 



if p^ + 1^"^ <l andO < p <l. 

In the next theorem we compute the equivalent of the variance of the number 
of specular points, under certain hypotheses on the random process and with 
the longuet-Higgins asymptotic. This result is new and useful for estimation 
purposes since it implies that, as fc ^ 0, the coefficient of variation of the 
random variable S tends to zero at a known speed. Moreover, it will also appear 
in a natural way when normalizing S to obtain a Central Limit Theorem. 

Theorem 7 Assume that the centered Gaussian stationary process W — {W{x) 
X G M} is 6— dependent, that is, T(z) = if \z\ > S, and that it has C^-paths. 
Then, as k ^ we have: 

Var(S') = 6*^ + 0(1). (22) 

k 

where 

/J , /2A7 25x4 



V^/2 V TT Vtt^As 
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J 



-s V2^(A2+r"(z)) 



dz. 



(23) 



the functions H and cf'^(z) have already been defined above, and 



a^z) 



T"'{zYT"{z) 



Al-r"2(z) 



Remarks on the statement. 



• The assumption that the paths of the process are of class imply that 
As < oo. This is well-known for Gaussian stationary processes (see for 
example [8,). 

• Notice that since the process is (5-dependcnt, it is also (5'-dcpcndcnt for any 
5' > 6. It is easy to verify that when computing with such a S' instead of 
S one gets the same value for 9. 

• One can replace the (S-dependence by some weaker mixing condition, such 
as 

|r(')(z)| < (const)(l + |z|)-" (0<i<4) 
for some a > 1, in which case the value of 6 should be replaced by: 

'•+°° ra2(z)iJ(p(z);0,0) 



2A4 

TT 



1 



1 A4 



2v/A2+r"(z) 



TT \/A2 



dz. 



The proof of this extension can be performed following the same lines as 
the one we give below, with some additional computations. 



Proof of the Theorem: We use the notations and computations preceding 
the statement of the theorem. 

Divide the integral on the right-hand side of (fT4)) into two parts, according as 
|x — yl > S or \x — y\ < 6, i.e. 

E(5(5-l))= / / •■•+// ... = /i+/2. (24) 

J J\x~v\>5 J J\x-y\<S 

In the first term, the (5— dependence of the process implies that one can 
factorize the conditional expectation and the density in the integrand. Taking 
into account that for each x e M, the random variables W"{x) and W'{x) are 
independent, we obtain for Ii: 

11 ^\W''{x)-k\)¥.{\W"{y)^k\)pw,^,-^{kx)pw'iy){ky)dxdy. 

J J\x-v\>S 

On the other hand, we know that W'{x) (resp. W"{x)) is centered normal with 
variance A2 (resp. A4). Hence: 

/i = [G(fc, V A4)J / / TT^'^^V -7; ^ dxdy, 

J J\x-y\>S II A2 J 
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To compute the integral on the right-hand side, notice that the integral over 
the whole x^y plane is equal to so that it suffices to compute the integral 
over the set \x — y\ < 5. Changing variables, this last one is equal to 



+00 



dx 



1 



■ exp 



1 



27rfc2 
5 



x-s 27rA2 

e^2" du 



1 k^{x^+y^) 



A2 



e ^ dv 



dy 



0(1), 



where the last term is bounded if k is bounded (in fact, remember that we are 
considering an approximation in which fc w 0). So, we can conclude that: 



1 



\x-y\>5 



2ttX 



■ exp 



ie{x^ + y^) 



A, 



dxdy 



1 



k\f\^ 



0[\) 



Replacing in the formula for 1\ and performing a Taylor expansion, we get: 



Let us now turn to I1. 



2A4 



1 



kJXrTi 



0{1) 



(25) 



Using Lemma [Hand the equivalences ([TSl) and pT|) . whenever \z\ = \x — y\ < 
S, the integrand on the right-hand side of (|14p is bounded by 

(const) [H{p{z);0, 0) + k^{x'^ + y^)] . 

We divide the integral I2 into two parts: 

First, on the set {{x,y) : \x\ < 26,\x — y\ < 6} the integral is clearly bounded 
by some constant. 

Second, we consider the integral on the set {{x,y) : x > 25,\x — y\ < 5}. 
(The symmetric case, replacing x > 2d by x < ~2S is similar, that is the reason 
for the factor 2 in what follows). 

We have (recall that z — x — y): 



/2 = 0(l) + 2/ / a^iz) H{p{z)■0,Q)+R2{p{z)■^,,v) 

l\x-y\<S,x>2S 



1 



2^^Xi-T"^{z) 



exp 



1 fc2(A2x2 + 2r"(a; - y)xy + \2y^) 



\l-T"^{x-y) 



dxdy 
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which can be rewritten as: 
r 



1 



V2^(A2+r"(z)) 

+ OC- ^ 



exp 



1 yz 



2^2 



A, 



2A2-r"(z)^A2 + r"(z) 



: exp 



(x- z/2)2 



/25 ^27r(A2 - r"(z)) 
In the inner integral we perform the change of variables 

V2k{x - z/2) 



A2-r"(z)) 



dx 



- r"(z) 



so that it becomes: 

1 /•+°° 1 



= exp ( T Idr 

27r 2 ^ 



1 1 



fcV2^ro 

where tq = 2V2fc(2(5 - z/2)/^A2 - r"(z) 



2^/2 A: 



0(1) 



dz 



(26) 



Notice that 0(1) in is uniformly bounded, independently of fc and z, since 
the hypotheses on the process imply that A2 — ^"{z) is bounded below by a 
positive number, for all z. 



We can now replace in the expression for I2 and we obtain 

/2=0(1)+ 



(27) 



To finish, put together ^ with (jM)) . (fT^ and ([T^ . 



Corollary 1 Under the conditions of Theorem^ as k ^ 0; 

The proof follows immediately from the Theorem and the value of the ex- 
pectation. 

The computations made in this section are in close relation with the two 
results of Theorem 4 in Kratz and Leon [TJ] . In this paper the random variable 
SP2{I) is expanded in the Wiener-Hermite Chaos. The aforementioned expan- 
sion yields the same formula for the expectation and allows obtaining also a 
formula for the variance. However, this expansion is difficult to manipulate in 
order to get the result of Theorem [71 



Let us now turn to the Central Limit Theorem. 
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Theorem 8 Assume that the process W satisfies the hypotheses of Theorem^ 
In addition, we will assume that the fourth moment of the number of approxi- 
mate specular points on an interval having length equal to 1 is hounded uniformly 
in k, that is 

E([^P2([0,1])]^) < (const) (28) 



Then, as k 0, 

S — /2A4 1 

' ' " ^ N{0,1), 



y/0/k 

where =^ denotes convergence in distribution. 
Remark. 



One can give conditions for the added hypothesis ([28| to hold true, which 
require some additional regularity for the process. Even though they are not 
nice, they are not costly from the point of view of physical models. For example, 
either one of the following conditions imply ([28]): 

• The paths x ~-+ W{x) are of class C^^. (Use Theorem 3.6 of [¥ with rn ~ 4, 
applied to the random process {W^'(x) : a; e K}. See also [IB]). 

• The paths x ^ W{x) are of class and the support of the spectral 
measure has an accumulation point: apply Exercice 3.4 of [4] to get the 
non-degeneracy condition. Proposition 5.10 of [4j and Rice formula (The- 
orem [2]) to get that the fourth moment of the number of zeros of W" (x) 
is bounded. 

Proof of the Theorem. Let a and f3 be real numbers satisfying the conditions 
1/2 < a < 1, a + P > 1, 2a + f3 < 2. It suffices to prove the convergence as k 
takes values on a sequence of positive numbers tending to 0. To keep in mind 
that the parameter is k, we use the notation 

S{k) ■.^S = SP2{M.) 

Choose k small enough, so that k^" > 2 and define the sets of disjoint 
intervals, for j = 0, ±1, . . . , ±[fc~'^]: 



= ((j - i)r"]5 + 5/2, j[fc-"]<5 - 5/2), 

=[j[k-'^]6-5/2,j[k-'^]S + S/2]. 
[.] denotes integer part. 

Notice that each interval has length [k~°']6 — S and that two neighboring 
intervals Uj are separated by an interval of length 5. So, the (5-dependence of 
the process implies that the random variables SP2{Uj), j = 0,±1, . . . ,±[k~^] 
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are independent. A similar argument applies to SP2{Ij), j = 0, ±1, . . . , ±[fc 
We denote: 

lil<[fc-"] 

Denote 

= (Var(^(fc)))"'/'«yfc7^ 
where the equivalence is due to Theorem [71 

We give the proof in two steps, which easily imply the statement. In the 
first one, we prove that 

Vk[S{k)-T{k)] 

tends to in the of the underlying probability space. 
In the second step we prove that 

VkT{k) 

is asymptotically standard normal. 

Step 1. We prove first that Vk[S{k) - T{k)] tends to in L^. Since it is 
non-negative, it suffices to show that its expectation tends to zero. We have: 

S{k) - T{k) = SP2{lf) + Z1 + Z2 

\j\<[k-i^] 

where 

Zi =5P2(-oo,-[fc-''].[fc-"]^ + (5/2), 
Z2 = SP2{[k-P].[k-'']S-5/2,+oo)). 
Using the fact that E(5P2 (-^)) ^ {const) Jj (p{kx / ^/X2)dx , one can show that 

VkE{S{k)~T{k)) < {const)k^/^\Yv{^^^^T^)+ C (p{kx/^)dx . 

which tends to zero as a consequence of the choice of a and (3. 

It suffices to prove that Var — T(fc)) ^ as fc ^ 0. Using indepen- 
dence: 

Var - T{k)) = ^ Var {SP2il!^)) + Var(Zi) + Var(Z2) 

b\<ik-n 

< Y HsP2{i!;){sP2{i^)~i)) (29) 

\3\<ik-n 

+ E{Zi{Zi - 1)) + E{Z2{Z2 - 1)) + E{S{k) ~ T{k)) . 
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We already know that Vj^ E{S{k) - T{k)) 0. Using the hypotheses of the 
theorem, since each 1^ can be covered by a fixed number of intervals of size one, 
we know that E(^SP2{Ij){SP2{Ij) — 1)) is bounded by a constant which does 
not depend on k and j. We can write 

J2 E{SP2iI^){SP2il!^)-l)) <{const)k^-^ 

\3\<lk-0] 

which tends to zero because of the choice of /3. The remaining two terms can 
be bounded by calculations similar to those of the proof of Theorem [T] 

Step 2. T{k) is a sum of independent but not equi-distributed random 
variables. To prove it satisfies a Central Limit Theorem, we use a Lyapunov 
condition based of fourth moments. Set: 

M™ e{ [SP2{U!f) - E{SP2{U!f 
For the Lyapunov condition it suffices to verify that 

Y,-^ Mj~>Q as fc^O, (30) 

b1<[fc-^] 

where 

\o\<[k-^\ 

To prove ([30]) . let us partition each interval into p — [fc^"] — 1 intervals 
Ii, ...Ip of equal size S. We have 

E{SPi + --- + SPpf = HSP,,SP,,SP,,SP,,), (31) 

1<Z1 ,22 ,23 ,'t4<p 

where SP^ stands for SP2{h) - E{SP2{h)) Since the size of all the intervals 
is equal to S and given the finiteness of fourth moments in the hypothesis, it 
follows that E{SP,,SP^,_SP^^SP^^) is bounded. 

On the other hand, notice that the number of terms which do not vanish in 
the sum of the right-hand side of pTjl is 0{p'^). In fact, if one of the indices in 
(*i, *2, *3, *4) differs more than 1 from all the other, then E^SPi-^SPi^SPi^SPi^^ 
vanishes. Hence, 



E 



SP2{U^) - E{SP2{U^))V < {const)k-^" 



so that E|j|<[/c-f] Mf = ©(fc-^^fc-zS). The inequality 2a + (3 < 2 implies Lya 
punov condition. | 
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3.3 Number of specular points without approximation 

We turn now to the computation of the expectation of the number of specular 
points SPi{I) defined by (P). This number of specular points is equal to the 
number of zeros of the process 

Z{x) := W'{x) - miix, W{x)) = 0, 

where 



x"^ - (hi - w){h2 -w) + ^/[x2 + {hi - u;)2][a;2 + {h2 - wY] 

■mi(X,W) — T- ; — . 

^ ' ' x{hi + h2-2w) 

Assume that the process : x G M} is Gaussian, centered, stationary, with 

Ao = 1. The process Z{t) is not Gaussian and we must use a generalization of 
Theorem [U namely Theorem 3.2 of |4] to get 

/'+00 

E{SPi{[a,b])) = dx E{\Z\x)\\Z{x) ^0,W{x) ^w) 

J a J ~oo 

1 „2 1 _^rjq(^^ 

e~~ =e dy^_ (32) 



For the conditional expectation in ((32)) . notice that 



Z\x) = W"{x) - W{x)) - ^{x, W{x))W'{x), 

so that under the condition, 

Z'{x) — W" {x)—K{x,w), where K{x,w) — — — {x,w))+— — {x,w))mi{x,w). 

ox ow 

Using that for each x, W"{x) and W'{x) are independent random variables 
and performing a Gaussian regression of W"{x) on W{x), we can write ([32]) in 
the form: 



E{SPi{[a,b])) 

pb r-\-oo 1 1 2 ( \ 

= dx E(|C-A2«;-i^(x,u;)|)-4=exp(-i(^«2 + ^^^i^))d^^ 

(33) 

where C is centered Gaussian with variance A4 — A|. Formula (j33p can still be 
rewritten as: 

E(5Fi([a,6])) 

= ^i^i '^L GKl)exp(-2(-^ + ^\^))^-, (34) 
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where 

, , X2W + K{x,w) 
m = in{x, w) = ^ — . 

Notice that in the integral is convergent as a ^ —00, b ^ +00 and 

that this formula is well-adapted to numerical approximation. 



3.4 Number of twinkles 

We give a proof of a result stated in [M] (part III pages 852-853). 

We consider Y{x,t) defined by ^ and we limit ourselves to the case in which 
W{x,t) is centered and stationary. If Y satisfies the conditions of Theorem |3l 
by stationarity we get 

E(r>V(/,T)) =T J^E{\dctY'{x,t)\\Y{x,t) = 0)pY(x,t)(0)dx. (35) 

Since Wxx and Wx are independent with respective variances 

^ + 00 /'+C30 



A40 — 



efi{d^,dT) X2o= / eKd^,dT), 

-00 J —00 



where jj. is the spectral measure of the stationary random field W{x,t). The 
density in ([35l) satisfies 



PY(.,t)(0) = (A2o)-^/V(A:a;(A2o)-^/2)(A4o)-^/V(fc(A4o)-'/') 
On the other hand 



Y'{x,t) = 



Wxx{x,t)-k Wxt{x,t) 
Wxxx{x,t) Wxxt{x,t) 



Under the condition Y{x,t) = 0, one has 

|det(Y'(x,t))| = \Wxt{x,t)Wxx.(x,t)l 

Computing the regression it turns out that the conditional distribution of the 
pair {Wxtix, t), Wxxxix, t)) under the same condition, is the one of two indepen- 
dent centered gaussian random variables, with the following parameters: 

A 1 A^ 
expectation - — k and variance A22 — for the first coordinate (36) 
A40 A40 

A4 A^ 
expectation - — kx and variance Ago — -r^' second coordinate (37) 

A20 A20 

It follows that: 



E 



(1 det(Y'(x, t))| I Y(a;, t) - o) - ^(^^' \/^22-^) •^(^fcx, ^A 



•^40 

60 — T 

A20 ' 
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Summing up: 

^E(r>V(R, T)) 




•^40 ^ k \/A2oA4o Afi + A/ 



A20 

(38) 



setting Xq :— Ago — This result is equivalent to formula (4.7) of [2] (part 
III page 853). 

3.5 Specular points in two dimensions 

We consider at fixed time a random surface depending on two space variables x 
and y. The source of light is placed at (0, 0, hi) and the observer is at (0, 0, ft.2)- 
The point (x, y) is a specular point if the normal vector n{x, y) — {—Wx, —Wy, 1) 
to the surface at (x, y) satisfies the following two conditions: 

• the angles with the incident ray / = (— —y, hi — W) and the reflected 
ray R = {—x, —y, /i2 — W) are equal (for short the argument (x, y) has 
been removed), 

• it belongs to the plane generated by / and R. 



Setting ai = hi — W and ri — x"^ + y"^ + ai, i = 1, 2, as in the one-parameter 
case we have: 

X Q2ri - air2 
= 2 I 2 ' 

y ayVi — airy , , 

When hi and /12 are large, the system above can be approximated by 

Wx — kx 

Wy = ky, (40) 
under the same conditions as in dimension 1. 

Next, we compute the expectation of SP2{Q), the number of approximate 
specular points in the sense of that are in a domain Q. In the remaining of 
this paragraph we limit our attention to this approximation and to the case in 
which {W{x,y) : {x,y) E R^} is a centered Gaussian stationary random field. 
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Let us define: 

Under very general conditions, for example on the spectral measure of y) : 

x,y gM.} the random field {Y{x, y) : x,y gM.} satisfies the conditions of Theo- 
rem [31 and we can write: 

E{SP2{Q))= f E(|detY'(a:,y)|)pY(..,)(0) dxdy, (42) 
Jq 

since for fixed (x, y) the random matrix Y'(x, y) and the random vector Y(x^ y) 
are independent, so that the condition in the conditional expectation can be 
erased. 



The density in the right hand side of has the expression 

exp -— — ^ 72-^(Ao2a;^ - 2Aiia;j/ + A2oj/^) 



P^ix,y}iO) ^P{w^,Wy){kx,ky) 
1 1 



2"' \/ A20A02 — 



2(A2oAo2 — 



(43) 

To compute the expectation of the absolute value of the determinant in the right 
hand side of (j42|) . which does not depend on x, y, we use the method of [6]. Set 
A := det Y'(a:, y) = (W,, - k){Wyy - k) - W^y. 



We have 



E(|A|) =E 



1 - cos(At) 



dt 



Define 
Then 



h{t) :- E [exp (it[{W.^ - k){Wyy - k) - W^y])] 
E(|A|) = ^(r-l-^^['^(^)]^ 



To compute h{t) we define 



A 



1/2 
1/2 
0-1 



and S the variance matrix of Wxx, Wyy, Wx 



vy-i 

A40 A22 A31 

A22 Ao4 A13 

A31 Ai3 A22 



(44) 



(45) 
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LetEV2AI]i/2 = P dwg(Ai, A2, A3)P'^ where P is orthogonal. Then by a 
diagonahzation argument 

h{t) = 6**'=' 

E(^exp[it{{AiZf-k{sn+S2i)Zi)+{A2Zi-k{su+S22)Z2)H^3Zi~k{sis+S23)Z3 



where {Zi, Z2, Z3) is standard normal and stj are the entries of E^/^P^. 

One can check that if ^ is a standard normal variable and r, are real 
constants, r > 0: 



— 2t / T/i 



where ^ 

= - arctan(2T), < < 7r/4. 

Replacing in pS)) . we obtain for 9^e[/i(t)] the formula: 
where, for j — 1, 2, 3: 



2 



1 + 4r2 ' 



mtm] = [n / 1 cos + (47) 

^,=l Jl + 4A2t2J ^ 



• fc) = exp 



2 1 + 4A2t2 



• fjit) — — arctan(2Aji), < ipj < 7r/4, 

' 3 l + 4A2t2 

Introducing these expressions in (j45p and using (j43p we obtain a new formula 
which has the form of a rather complicated integral. However, it is well adapted 
to numerical evaluation. 

On the other hand, this formula allows us to compute the equivalent as 
fc ^ of the expectation of the total number of specular points under the 
longuet-Higgins approximation. In fact, a first order expansion of the terms in 
the integrand gives a somewhat more accurate result, that we state as a theorem: 



Theorem 9 

E(5P2(M')) = ^ + 0(l) (48) 
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where 



m2 = p dt 

+- 1 - 2-y^[UU {^jV^+A])] (1 ~ B,B, - ij2^3 - B^B,) 

(49) 

where 







Notice that m2 only depends on the eigenvalues Ai , A2 , A3 and is easily 
computed numerically. 

In Flores and Leon 10 a different approach was followed in search of a for- 
mula for the expectation of the number of specular points in the two-dimensional 
case, but their result is only suitable for Montecarlo approximation. 

We now consider the variance of the total number of specular points in 
two dimensions, looking for analogous results to the one-dimensional case (i.e. 
Theorem [7] and its Corollary [1]) , in view of their interest for statistical appli- 
cations. It turns out that the computations become much more involved. The 
statements on variance and speed of convergence to zero of the coefficient of 
variation that we give below include only the order of the asymptotic behavior 
in the longuct-Higgins approximation, but not the constant. However, we still 
consider them to be useful. If one refines the computations one can give rough 
bounds on the generic constants in Theorem [TOl and Corollary [2] on the basis of 
additional hypotheses on the random field. 

We assume that the real- valued, centered, Gaussian stationary random field 
{W^(x) : X g M^} has paths of class C^, the distribution of W'{0) does not 
degenerate (that is Var(VF'(0)) is invertible). Moreover, let us consider W"{0), 
expressed in the reference system xOy of as the 2x2 symmetric centered 
Gaussian random matrix: 



VF"(0) 

The function 



z A(z) = det [ Var {W" {0)z)] , 

defined on z = (21,22)"^ € R^, is a non-negative homogeneous polynomial of 
degree 4 in the pair zi, Z2. We will assume the non-degeneracy condition: 

min{A(z) : ||z|| = 1} = A > 0. (50) 

Theorem 10 Let us assume that {M^(x) : x G R^} satisfies the above condtions 
and that it is also 6-dependent, S > 0, that is, K(W{x)W{y)) = whenever 



22 



I|x-y|l ><^. 

Then, for k small enough: 

Var(5P2(M^)) < ^, 

where L is a positive constant depending upon the law of the random field. 

A direct consequence of Theorems [5] and [TU] is the fohowing: 

Corollary 2 Under the same hypotheses of Theorem \1(A for k small enough, 
one has: 

^ , < Life 

E(5P2(R2)) - 

where L\ is a new positive constant. 

Proof of Theorem[T^ For short, let us denote T — SP2(M.'^). We have: 

Var(r) ^ E{T{T - 1)) + E(T) - [E{T)]^ (51) 

We have aheady computed the equivalents as fc ^ of the second and third 
term in the right-hand side of ([51]) . Our task in what follows is to consider the 
first term. 

The proof is performed along the same lines as the one of Theorem [3 but 
instead of applying Rice formula for the second factorial moment of the num- 
ber of crossings of a one-parameter random process, we need Theorem |4] for 
dimension d = 2. We write the factorial moment of order m = 2 in the form: 



E(T(T- 1)) 

= // E(|detY'(x)||detY'(y)|| Y(x) =0,Y(y) =0)pY(x).Y(y)(0,0) dxdy 



... dxdy + / / ••■ dxdy = Ji + J2. 

X— y||>(5 J J \\x—y\\<5 

For Ji we proceed as in the proof of Theorem [7j using the (5-dependence and 
the evaluations leading to the statement of Theorem [SI We obtain: 

T ^ 0(1) 

•^1-^ + ^- (52) 

Let us show that for small k, 

J2 = (53) 
In view of ([ST]) . ([48[) and ([52[l this suffices to prove the theorem. 
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We do not perform all detailed computations. The key point consists in 
evaluating the behavior of the integrand that appears in J2 near the diagonal 
X = y, where the density PY(x),Y(y)(0, 0) degenerates and the conditional ex- 
pectation tends to zero. 

For the density, using the invariance under translations of the law of VF'(x) : 
X e M^, we have: 

PY(x),Y(y)(0, 0) = Pw'(^).W'{y) (^^' ^y) 

= PW"(0),H"(y-x) (fcx, fcy) 

= PH"(o),[iV'(y-x)-H"(o)](fcx, k{y - x)). 

Perform the Taylor expansion, for small z = y — x G : 

P^'(z) = VK'(O) + M^"(0)z + 0(||z||2). 

Using the non-degeneracy assumption ([SO)) and the fact that VF'(O) and W"{Q) 
are independent, we can show that for x, z e M^, ||z|| < 5: 

PY(x).Y(y)(0,0) < ^exp [- Cafc^dlxll - C^f] 

where Ci, C2, C3 are positive constants. 

Let us consider the conditional expectation. For each pair x, y of different 
points in M^, denote by r the unit vector (y — x)/||y — x|| and n a unit vec- 
tor orthogonal to r. We denote respectively by (9tY, 9ttY, 9nY the first and 
second partial derivatives of the random field in the directions given by r and n. 

Under the condition 

Y(x) = 0,Y(y) -0 

we have the following simple bound on the determinant, based upon its definition 
and RoUe's Theorem applied to the segment [x, y] = {Ax + (1 — A)y}: 



|detY'(x)| < |15.Y(x)|l|15„Y(x)|l < |ly-x|l sup |ia..Y(s)|l |ia„Y(x)|l (54) 

se[x,y] 



So, 

E(|detY'(x)||detY'(y)|| Y(x) = 0,Y(y) = 



< l|y-x|pE 



sup ||a,,Y(s)||2||a„Y(x)||||a„Y(y)|| W'{^) = fcx,iy'(y) = fcy 

se[x,y] 

sup ||a,,Y(s)||2||a„Y(o)|||ia„Y(z)|i w'{q) ^ fcx, ^'(^) - ^'(Q) 

se[o,z] 
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where the last equahty is again a consequence of the stationarity of the random 
field {W^(x) : x e M.^}. 

At this point, we perform a Gaussian regression on the condition. For the 
condition, use again Taylor expansion, the non-degeneracy hypothesis and the 
independence of W(0) and W"{0). Then, use the finiteness of the moments of 
the supremum of bounded Gaussian processes (see for example [4], Ch. 2), take 
into account that |jz|| < S to get the inequality: 

E(|detY'(x)||detY'(y)|| Y(x) = 0, Y(y) = o) < C4 ||z||2(l + A;||x||)' (55) 

where C4 is a positive constant. Summing up, we have the following bound for 

J2<CiCinS^ f (1 + A:|lx|l)^exp[- 6*2^2(11x11 -C3)"]dx 

.+00 (56) 
= C1C4 2Tr^d^ / (1 + kpyexp[^C2k^{p-C3f]pdp 
Jo 

Performing the change of variables w = kp, (|53p follows. | 



4 The distribution of the normal to the level 
curve 

Let us consider a modeling of the sea W{x, y, t) as a function of two space vari- 
ables and one time variable. Usual models are centered Gaussian stationary 
with a particular form of the spectral measure /x that we discuss briefly below. 
We denote the covariance by r(a;, y, t) = E(VF(0, 0, Q)W{x, y,t)). 

In practice, one is frequently confronted with the following situation: several 
pictures of the sea on time over an interval [0, T] are stocked and some properties 
or magnitudes are observed. If the time T and the number of pictures are large, 
and if the process is ergodic in time, the frequency of pictures that satisfy a 
certain property will converge to the probability of this property to happen at 
a fixed time. 

Let us illustrate this with the angle of the normal to the level curve at a 
point "chosen at random" . We consider first the number of crossings of a level 
u by the process W{-,y,t) for fixed t and y, defined as 

= #{x : < X < Mi; W{x,y,t) = u}. 

We are interested in computing the total number of crossings per unit time when 
integrating over y g [0, M2] i.e. 

i dtj^ 7V[^l^f)(u)dy. (57) 
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If the ergodicity assumption in time holds true, we can conclude that a.s. 



1 f^'^ Ar'^i-.y.t)), ^, AT W( AT'^i-'O.O)), n\ M1M2 / A2OO -i- 



where 



A"A(^Aj(i/i(A2;, Aj,, Af) 



are the spectral moments. 



Hence, on the basis of the quantity (|57p for large T, one can make inference 
about the value of certain parameters of the law of the random field. In this 
example these are the spectral moments A200 and Aqoo- 

If two-dimensional level information is available, one can work differently be- 
cause there exists an interesting relationship with Rice formula for level curves 
that we explain in what follows. 

We can write (x = (x, y)): 

W'{yL,t) = ||M/'(x,i)||(cose(x,i),sine(x,i))^. 
Instead of using Theorem [Tl we can use Theorem [Bl to write 



Ah p 

KAitf^ dy] - E [ / I cos e(x, 0) I da 

' ' ' ■JCQ(a,u) 



CT2(Q) /A200 /cox 
= \ - e , (58) 

Ti" V Aooo 

where Q = [0, Mi] x [0, M2]. We have a similar formula when we consider sec- 
tions of the set [0, Mi] x [0, M2] in the other direction. In fact ([55)1 can be 
generalized to obtain the Palm distribution of the angle 8. 

Set hg^_0r^ — e^Jj and for — tt < 61 < 62 < define 
F(^2) - F{ei) : = e((7i({x e Q : M^(x,0) = u ; 0i < e(x,s) < ^2}) 



= E / V,e.(e(x,s))rfai(x)ds (59) 

Denoting A = A200A020 — Ano and assuming a2{Q) = 1 for ease of notation, we 
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readily obtain 



F{92)-F{e^) 



2 



(2^)3/2 (A) 1/2 VA^y„. 



2 



e 




p2 exp( 




where A_ < A+ are the eigenvalues of the covariance matrix of the random vector 
{dxW{0,0,0),dyW{0,0,0)) and k is the angle of the eigenvector associated to 
7+. Remarking that the exponent in the integrand can be written as 
1/A_(1 — 7^szn^((^ — k)) with 7^ 1 — A+/A_ and that 



From this relation we get the density g{ip) of the Palm distribution, simply by 
dividing by the total mass: 



Here /C is the complete elliptic integral of the first kind. This density character- 
izes the distribution of the angle of the normal at a point chosen "at random" 
on the level curve. 

In the case of a random field which is isotropic in {x,y), we have A200 = -^020 
and moreover Ano = 0, so that g turns out to be the uniform density over 
the circle (Longuet-Higgins says that over the contour the "distribution" of the 
angle is uniform (cf. (TS], pp. 348)). 

Let now W = {W{x,t) : t e K+,x = (x, y) G K^} be a stationary zero 
mean Gaussian random field modeling the height of the sea waves. It has the 
following spectral representation: 




it is easy to get that 





(60) 
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where A is the manifold {Af + = u;'^} (assuming that the acceleration of 
gravity g is equal to 1) and M is a random Gaussian orthogonal measure defined 
on A (see [13] )■ This leads to the following representation for the covariance 
function 



Tix,y,t)= / e*(^i-+^^^+-*)/(Ai,A2,c^)a2(dF) 

J A 

r2Tr 

^i{uj^x cos ^+uj^y sin ip+ujt) / 



where, in the second equation, we made the change of variable Ai — uj^ cos (p, 
A2 — oj^sin^s and G{ip,ui) = f{uj^ cos. ip^uP' sin ip,uj)2bj'^ . The function G is 
called the "directional spectral function". If G does not depend of Lp the ran- 
dom field W is isotropic in x, y. 

Let us turn to ergodicity. For a given subset Q of and each i, let us define 

At = (j{W{x, y,t) -.T >t;{x,y) e Q} 

and consider the cr-algebra of i-invariant events A ~ Pl-^*- We assume that 
for each pair {x,y), T(x,y,t) — s- as i — s- +00. It is well-known that under 
this condition, the ci-algebra A is trivial, that is, it only contains events having 
probability zero or one (see for example [8], Ch. 7). 

This has the following important consequence in our context. Assume further 
that the set Q has a smooth boundary and for simplicity, unit Lebesgue measure. 
Let us consider 

Z{t) = / H{x,t)dai{x.), 

JcQ(u,t) 

where iJ(x,i) = 7^(VF(x, i), VM/(x, t)) , where VW = {W^,Wy) denotes gradi- 
ent in the space variables and Ti. is some measurable function such that the inte- 
gral is well-defined. This is exactly our case in ([59)1 . The process {Z{t) : t &R}is 
strictly stationary, and in our case has a finite mean and is Riemann-integrable. 
By the Birkhoff-Khintchinc crgodic theorem (|F page 151), a.s. as T ^ +00, 



l l\{s)ds^EB[ZiO)], 



where B is the cr-algebra of ^-invariant events associated to the process Z{t). 
Since for each t, Z{t) is -measurable, it follows that i3 C -4, so that Eb[Z(0)] = 
£[^(0)]. On the other hand. Rice's formula yields (take into account that sta- 
tionarity of W implies that W^(0,0) and VVF(0, 0) are independent): 

E[Z(0)] = E[H{u, Viy (0, 0)) 1 1 VT4^(0, 0)\\]pwio,o) (u). 
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We consider now the CLT. Let us define 

ft 



Z{t)^\j^ [Z(s)-E(Z(0))]ds, 



In order to compute second moments, we use Rice formula for integrals over 
level sets (cf. Theorem [6]), applied to the vector- valued random field 

X(xi,X2,Si,S2) = (VF(xi,Si),iy(x2,S2))^. 

The level set can be written as: 

Cq2(u, u) = {(xi, X2) e QxQ : X(xi, X2, si, S2) = (m, u)} for < si < i, < S2 < t. 
So, we get 



Var Z(t) = - / (1 - s)ds, 



where 

/(m,s)= / E[i?(xi,0)iI(x2,s)|lVW^(xi,0)|l|lVVK(x2,s)|l VF(xi, 0) = m ; VK(x2, s) = « 

X Pw{^^,o),w{^2,s){u,u)dxldx2 - VW^(0,0))||VVK(0,0)||]pvy(o,o)(M))' 

Assuming that the given random field is time-5-dependent, that is, 
T{x,y,t) = V {x,y), whenever t > 5, we readily obtain 

5 

2/ 



iVarZ(t)-^2/ I{u, s)ds := a'' (u) as t ^ 00. 
"'0 

Using now a variant of the Hoeffding-Robbins Theorem [TT] for sums of 5- 
dependent random variables, we get the CLT: 

^ftZ{t) ^ A^(0,cr2(u)). 



5 Numerical computations 

Validity of the approximation for the number of specular 
points 

In the particular case of stationary processes we have compared the exact ex- 
pectation given by ([32|l with the approximation pO|l . 

In full generality the result depends on hi,h2, A4 and A2 . After scaling, we 
can assume for example that A2 = 1. 

The main result is that, when hi « /i2, the approximation (|10[) is very sharp. 
For example with the value (100, 100, 3) for {hi, h2, A4), the expectation of the 
total number of specular points over M is 138.2; using the approximation (llip 
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the result with the exact formula is around 2.10~^ larger but it is almost hidden 
by the precision of the computation of the integral. 

If we consider the case (90,110,3), the results are respectively 136.81 and 
137.7. 

In the case (100, 300, 3), the results differ significantly and Figure [1] displays 
the densities ^ and PH)) 




Figure 1: Intensity of specular points in the case hi = 100, /12 — 300, A4 = 3. 
In solid line exact formula, in dashed line approximation (jlOp 



Effect of anisotropy on the distribution of the angle of the 
normal to the curve 

We show the values of the density given by (pO)) in the case of anisotropic 
processes 7 = 0.5 and k = it/ A. Figure [2] displays the densities of the Palm dis- 
tribution of the angle showing a large departure from the uniform distribution. 



Specular points in dimension 2 

We use a standard sea model with a Jonswap spectrum and spread function 
cos(20). It corresponds to the default parameters of the Jonswap function of 
the toolbox WAFO [TH]. The variance matrix of the gradient is equal to 

10- (^^^ M 

\ Q 81 y 
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Figure 2: Density of the Palm distribution of the angle of the normal to the 
level curve in the case 7 = 0.5 and n — 7r/4 



and the matrix E of Section 13.51 is 

/ 9 3 

E = 10"-* 3 11 

\ 3 

The spectrum is presented in Figure [3] 




Figure 3: Directional Jonswap spectrum as obtained using the default options 
of Wafo 
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The integrand in (|^^ is displayed in Figure [3] as a function of the two space 
variables x, y. The value of the asymptotic parameter m2 defining the expansion 
on the expectation of the numbers of specular points, see 1^5)) . is 2.52710"^. 




Figure 4: Intensity function of the specular points for the Jonswap spectrum 

The Matlab programs used for these computations are available at 
\protect\vrule widthOpt\protect\href {http : //www. math. univ-toulouse . f r/\string~azais/prog/pri 

6 Application to dislocations of wavefronts 

In this section we follow the article by Berry and Dennis [5] . As these authors, we 
are interested in dislocations of wavefronts. These are lines in space or points in 
the plane where the phase x, of the complex scalar wave ip{:ii, t) = /3(x, t)e^^^^'*\ 
is undefined, (x — (xi,X2)) is a two dimensional space variable). With respect 
to light they are lines of darkness; with respect to sound, threads of silence. 
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It will be convenient to express ip by means of its real and imaginary parts: 

fp{x,t) = ^(x,i) +ir/(x,t). 

Thus the dislocations are the intersection of the two surfaces 

C(x,t) = ??(x,t) = 0. 

We assume an isotropic Gaussian model. This means that we will consider the 
wavefront as an isotropic Gaussian field 

V'(x, t)= f exp (z[(k . x) - c\m){^^)'^'dW{k), 

where, k = (fci,fc2), |k| = \/k1 + n(fc) is the isotropic spectral density 
and W = {Wi + iW2) is a standard complex orthogonal Gaussian measure 
on M^, with unit variance. Here we are interested only in t = and we put 
^(x) :=C(x,0) and ?7(x) :=77(x,0). 

We have, setting fc = |k| 

€(x) = / cos((k.x))(ffl)V2dW^,(k)- / sin((k.x))(M)i/2rfiy2(k) (gi) 

and 

r?(x)=/ cos((k-x))(ffl)i/2dW'2(k)+ / sin((k-x))(ffl)i/2rf^^(k) (52) 
The covariances are 

/■OO 

E [C(x)C(x')] = E [7?(x)77(x')] = p(|x - x'l) := / Jo(fc|x - x'|)n(fc)dfc (63) 

where Ju{x) is the Bessel function of the first kind of order u. Moreover 
E[e(ri)r7(r2)] =0. 

Three dimensional model 

In the case of a three dim ensional Gau ssian field, we have x = {x\,X2,Xi), 
k = {ki,k2, kz),k = |k| = ■\/fcfT"fcf+l!f and 

V(x) = / exp (i[(k • x)]) (IHP) V2rfw^(k). 

In this case, we write the covariances in the form: 

E[^(r0^(r2)] = '-^^^l^U{k)dk. (64) 

The same formula holds true for the process r] and also E [^(ri)77(r2)] = for 
any ri,r2, showing that the two coordinates are independent Gaussian fields . 
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6.1 Mean length of dislocation curves, mean number of 
dislocation points 

Dimension 2: Let us denote {Z(x) : x e M^} a random field with values in R^, 
with coordinates ^(x),?7(x), which are two independent Gaussian stationary 
isotropic random fields with the same distribution. We are interested in the 
expectation of the number of dislocation points 

d2 :=E[#{xe5:e(x)=77(x)=0}], 

where 5* is a subset of the parameter space having area equal to 1. 

Without loss of generality we may assume that Var(^(x)) Var(?7(x)) = 1 
and for the derivatives we set A2 = Var(?7i(x)) ~ Var(^i(x)), i — 1,2. Then, 
using stationarity and the Rice formula (Theorem [3]) we get 

d2 = E[| det(Z'(x))|/Z(x) = 0]pz(x)(0), 

The stationarity implies independence between Z(x) and Z'(x) so that the con- 
ditional expectation above is in fact an ordinary expectation. The entries of 
Z'(x) are four independent centered Gaussian variables with variance A2, so 
that, up to a factor, | det(Z'(x))| is the area of the parallellogram generated by 
two independent standard Gaussian variables in R^. One can easily show that 
the distribution of this volume is the product of independent square roots of 
a x^(2) and a x^(l) distributed random variables. An elementary calculation 
gives then: E[| det(Z'(x))|] = A2. Finally, we get 



This quantity is equal to in Berry and Dennis 6] notations, giving their 
formula (4.6). 

Dimension 3: In the case, our aim is to compute 

ds = E[/:{x e S : ^(x) = ?7(x) = 0}] 

where is a subset of M.^ having volume equal to 1 and £ is the length of the 
curve. Note that is denoted by d [6j. We use the same notations and remarks 
except that the form of the Rice's formula is (cf. Theorem O 

d3 = :^E[(detZ'(x)Z'(xf 

Again 

E[(det(Z'(x)Z'(x)^)i/2] = A2E(P^), 

where V is the surface area of the parallelogram generated by two standard 
Gaussian variables in . A similar method to compute the expectation of this 
random area gives: 



E{V) = E( Vx^) X EiVxm) = 
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Leading eventually to 

= — . 

TT 

In Berry and Dennis' notations 6_ this last quantity is denoted by ^ giving 
their formula (4.5). 

6.2 Variance 

In this section we limit ourselves to dimension 2. Let S be again a measurable 
subset of having Lebesgue measure equal to 1. The computation of the 
variance of the number of dislocations points is performed using Theorem 2] to 
express 

E(7V|(0)(7V|(0)-1)) - / A(s,,S2)dsids2. 

We assume that {Z(x) : x e M^} satisfies the hypotheses of Theorem |4] for 
m = 2. Then use 



Var (iV|(0)) = E(7V|(0)(7V|(0) - l)) + ^2 - rf: 



Taking into account that the law of the random field is invariant under transla- 
tions and orthogonal transformations of M^, we have 

A(si,S2) = A((0,0),(r,0)) =A{r) whith r = ||si - S2II , 

The Rice's function A{r)) , has two intuitive interpretations. First it can be 
viewed as 

A{r) = lim-ljE[iV(i?((0,0),e)) x iV(i?((r, 0), e))] . 

Second it is the density of the Palm distribution (a generalization Horizontal 
window conditioning of [8]) of the number of zeroes of Z per unit of surface, 
locally around the point (r, 0) given that there is a zero at (0, 0). 
A{r)/d2 is called "correlation function" in [6]. 

To compute A{r), we put fi,^2,'7i?72 for the partial derivatives of ^,77 with 
respect to first and second coordinate. 

and 

Air) = E[| det Z'(0, 0) det Z'(r, 0)||Z(0, 0) = Z(r, 0) = 02]pz(o,o),z(r,o)(04) 
= E[| (^772 - i2m) (0, 0) {^im - 6^1) (r, 0)) I |Z(0, 0) = Z(r, 0) = O2] 

PZ(0,0),Z(r,0)(04) (65) 

where Op denotes the null vector in dimension p. 
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The density is easy to compute 

1 



PZ(0,0),Z(r,0)(04) 



(27r)2(l-p2(^)) 



/•OO 

, where p{r) — i Jo{kr)Il{k)dk. 
Jo 



We use now the same device as above to compute the conditional expectation 
of the modulus of the product of determinants, that is we write: 



1 



w 



+ 00 



(1 - cos{wt)r^dt. 



(66) 



and also the same notations as in [6] 



C := 


= P{r) 


E = 


P'ir) 


H = 


-E/r 


F = 




Fo = 


-p"(o) 



The regression formulas imply that the conditional variance matrix of the vector 

W = (a (0) , 6 1'^, 0) , 6 (0) , 6 (r, 0) , 771 (0) , 7?! (r, 0) , 772 (0) , 772 (r, 0)' 
is given by 



with 



A 



S = Diag 



A,B,A,B 



F 



B 



E^C 



F- 



E^C 



Fo H 
H Fo 



Using formula ([66]) the expectation we have to compute is equal to 



1 /■^°° , 



+00 



i-^nii,o)-^r(-ti,o)-ir(o,<2)-^T(o,-t2) 



+ \Tit,,t2) + ^T(-ii, t2) + ^T{h, -h) + ir(-ti, -h)] (67) 



where 
with 



T{ti,t2) = E[exp {i{witi + W2t2)) 



wi = Ci{0)m{0) ~ r/i (0)6(0) = W1W7 - W3W5 
7^2 = 6(7-, 0)772 (r,0) - 771 (r, 0)6 (r,0) = W2W8 - W4 We. 
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T(<i, ia) = ]E( exp(iW'^HW)) where W has the distribution N{0, S) and 



n = 



V 



A standard diagonaUzation argument shows that 














V 




















-V 








V 















1 


■ h 









2 





^2 . 





r(ti,t2) =E(exp(iW^HW)) =E(exp(i^Aje|)), 

where the ^j's are independent with standard normal distribution and the 
are the eigenvalues of E^/^TiS^/^. Using the characteristic function of the 
distribution: 



E(exp(iW^HW)) = " 2iAj' 
i=i 



-1/2 



(68) 



Clearly 



and 












M 
















-M 




















with M = A^/^VB^/^. 

Let A be an eig envalue of E^/^^Si/^ It is easy to check that A is an eigen- 
value of AiAi^ . Respectively if and A| are the eigenvalues of A4A4^ , those 
of aj.g ±Ai (twice) and ±A2 (twice). 

Note that A? and A^ are the eigenvalues of MM'^ = A^/'^VBVA^^'^ or 
equivalently, of VBVA. Using 



E(exp(iW^7^W)) = (1+4(A2+a2)+16A2a 
where 



( 1 +4<r [VBVA) + 1 6 det (PSPyt) ) ^ 



VBVA = - 

4 



ii-Fo(^o ~ ) + tit2H{F — ) t\F(){F — ) + tit2H{Fo — 
tit2H{Fo — Y^T^-) + *2^o(^ ^ T^S^) tit2H{F — jt:^) + ^2-^0(^0 — T^J^-) 



So, 



4ir(2?S2?^) = (t2+t2)j^^(^^ 



^2 E'^C 
-^—^) + 2ht2H{F-^^) 



(69) 



1 -C2' 

16 det(I?SP^) = K - H'] [{Fo - -{F- ^%^f] (70) 
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giving 



T{h,t2) =E(cxp(zW^HW)) 



(l + (4 + tl)Fo{Fo - + 2hhH{F - - 



+ t 



C2' 



\ii [El {[E, - -{E- 



(71) 



Performing the change of variable t' — \J~K\t with A\ — Eo{Eq — j^fjj) the 
integral ((67|) becomes 



/• + 00 p + cc 

71"^ 7-00 J-oo 



2 



1 1 1 r 

^"i + t2i + <2+"2lTT(t!' 



1 



1 



4i 

r2 



+ tl) - 2^2^1*2 + tltlZ 1 + (i? + ii) + 2A2iit2 + tltlZ - 

+00 /' + CX2 



— / dti I a^2f'i ^2 



1 - 



l + (i?+t2)+t2<2^ 



1 + ^1 1 + ^2 + (<2 + ^2) + ^2^2^^ _ 4^2^2i2 



(72) 



where 



A 



H F(1-C^)-E^C 



7 - ^n- g 



2 - p^^ Fo(l-C')-E-^ 

In this form, and up to a sign change, this result is equivalent to Formula (4.43) 
of [6] (note that = F in [6]). 

In order to compute the integral ([7^ . first we obtain 



J -QO ''2 -"^ T~ ''2 



'-co ''2 ^ < ''2 

We split the other term into two integrals, thus we have for the first one 

1 1 



i rir 

2looti^l + (i? 



+ 4) - 2A2tit2 + t\tlZ 1 + tl 



dU 



1 



2(l + <2) 

1 



1 (I + tlZ)tl - 2A2ht2 
41 + 11- 2A2tit2 + (1 + tlZ)ti ' 

1 to — '^Z\t\t2 



rdU 



2(1 + t?) J-oo 4-2Z,ht2 



-dto 



where Z2 

Similarly for the second integral we get 
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1 f°° j_r 1 1_. 

2{l + tl)j_^tl 



tl 4 + 2Zitit2 + Z2 



dt 



2 



' ~ ~ 2(1 + i?) f _ 2Z^tit2 + Z2 4 + 2^1*1*2 + ^2 



J /-^ _ 

+ ^l) J-oo ^2 



dU 



(1 + t?) 4 + 2(^2 - 2Z2tf + z| 

1 ^(Z2 - 2Zltl) 

i^+tl) Z2^{Z2-Zlti)' 

In the third hne we have used the formula provided by the method of residues. 
In fact, if the polynomial X'^ — SX + P with P > has not root in [0, 00), then 



/. 



t bt /p(_5 + 2VP) 



In our case 7 = -{Z2 - 'iZftf), S = -2(^2 - 2Zftf) and P = Z|. 
Therefore we get 



4n^l~C^)J_^tl^' {l + ti)Z2J(Z^^ZM-y 



Acknowledgement 



This work has received financial support from European Marie Curie Network 
SEAMOCS. 



References 

[1] R. J. Adler, The Geometry of Random Fields, Wiley, (1981). 

[2] R. J. Adler and J. Taylor, Random Fields and Geometry. Springer, (2007). 

[3] J-M. Azai's, J. Leon and J. Ortega, Geometrical Characteristic of Gaussian 
sea Waves. Journal of Applied Probability , 42,1-19. (2005). 

[4] J-M. Azai's, and M. Wschebor, Level set and extrema of random processes 
and fields, Wiley (2009). 

[5] J-M. Azai's, and M. Wschebor, On the Distribution of the Maximum of a 
Gaussian Field with d Parameters, Annals of Applied Probability, 15 (lA), 
254-278, (2005). 



39 



[6] M.V. Berry, and M.R. Dennis, Phase singularities in isotropic random waves, 
Proc. R. Soc. Land, A, 456, 2059-2079 (2000). 

[7] E. Cabana, Esperanzas de Integrales sobre Conjuntos de Nivel aleato- 
rios. Actas del 2 Congreso Latinoamericano de Probabilidad y Estadistica 
Matematica, Editor: Sociedad Bernoulli seccion de Latinoamerica, Spanish 
, Caracas, 65-82 (1985). 

[8] H. Cramer and M.R. Leadbetter, Stationary and Related Stochastic Pro- 
cesses, Wiley (1967). 

[9] H. Federer, Geometric Measure, Springer (1969). 

[10] E. Floras and J.R. Leon, Random seas. Levels sets and applications, 
Preprint (2009). 

[11] W. HoefFding and H. Robbins, The Central Limit Theorem for dependent 
random variables, Duke Math. J. 15 , 773-780,(1948). 

[12] M. Kratz and J. R. Leon, Level curves crossings and applications for Gaus- 
sian models, Extremes, DOI 10.1007/sl0687-009-0090-x (2009). 

[13] P. Kree and C. Soize, Mecanique Aletoire, Dunod (1983). 

[14] M. S. Longuet-Higgins, Reflection and refraction at a random surface. I, II, 
III, Journal of the Optical Society of America, vol. 50, No. 9, 838-856 (1960). 

[15] M. S. Longuet-Higgins, The statistical geometry of random surfaces. Proc. 
Symp. Appl. Math., Vol. XIII, AMS Providence R.I., 105-143 (1962). 

[16] Nualart, D. and Wschebor, M., Integration par parties dans I'espace de 
Wiener et approximation du temps local, Prob. Th. Rel. Fields, 90, 83-109 
(1991). 

[17] S.O. Rice, (1944- 1945). Mathematical Analysis of Random Noise, Bell Sys- 
tem Tech. J., 23, 282-332; 24, 45-156 (1944-1945). 

[18] WAFO-group . WAFO - A Matlab Toolbox for Analysis of Random Waves 
and Loads - A Tutorial. Math. Stat., Center for Math. Sci., Lund Univ., Lund, 
Sweden. ISBN XXXX, URL|http://www.maths.lth.se/matstat/wafo.(2000[) 

[19] M. Wschebor, Surfaces Aleatoires. Lecture Notes Math. 1147, Springer, 
(1985). 

[20] U. Zahle, A general Rice formula. Palm measures, and horizontal-window 
conditioning for random fields, Stoc. Process and their applications, 17, 265- 
283 (1984). 



40 



